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Abstract. In this paper a new class of modified-Hessian equations, closely related to the Optimal Transportation 
Equation, will be introduced and studied. In particular, the existence of globally smooth, classical solutions of 
these equations satisfying the second boundary value problem will be proven. This proof follows a standard 
method of continuity argument, which subsequently requires various a priori estimates to be made on classical 
solutions. These estimates are modifications of and generalise the corresponding estimates for the Optimal 
Transportation Equation, presented in ITW06I . Of particular note is the fact that the global C 2 estimate 
contained in this paper makes no use of duality with regard to the original equation. 



1. Introduction 

In this paper the global regularity of classical solutions solving 



^y~' [D 2 u- D 2 x c(-,T U )} = B(-,u), in ft- (1.1a) 
associated with the second boundary value problem 

T (ft-) = ft+, (1.1b) 

for given c, B, T u and ft* c R" will be studied. Specifically, ft* are smooth domains in R" and 

i 

5 n \- 



5, 



[M] := d ni (Ai A n ), 

'S„(Ai A„)\^ 



S<(Ai Xn) 

with Ai, . . . , A„ being the eigenvalues of the n x n symmetric matrix M and Se denoting the £-th elementary 
symmetric function given by 



Si(Xi A n )= J2 [ti X J 

l<k<—<k<n \j=h 

for any integer I < n. The structure of T u and the precise conditions on ft*, c and B will be given in Section 

m 

The structure of (|l.la|) - (|l.lb|) is reminiscent of the Monge-Ampere type equations emerging from optimal 
transportation. Indeed, the conditions required and techniques used in this paper are reflective of those in 
|TW06| for the optimal transportation case, but all of the corresponding a priori estimates required to prove 
existence of globally smooth solutions in the current scenario need modification. In particular, (|l.la|) does not 
have a duality property like the optimal transportation equation, which is one of the obstacles preventing results 
from being carried over from optimal transportation to the current scenario. 

The author would like to thank Neil Trudinger and Xu-Jia Wang for their discussions on and advisement 
regarding the research contained in this paper. 
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2. Motivation and Main Results 

Before stating the main results of the paper, the mathematical motivation to study (|l.la|) - (|l.lb|) will be 
presented. Specifically, the class of boundary value problems depicted in Section Q] represent the only class of 
problems from which the estimates in |TW 06] can be adapted. 

2.1. The Structure of F. One can regard (|l.la|) as having the general form 

F [D 2 u - A(- , u, Du)] = 8(- , u), in n~, (2.1) 

where A is a given nx n matrix-valued function and 6 is a given scalar-valued function, defined on xlxl" 
and Q~ x R respectively. Subsequently, the left-hand side of (|2.ip can be represented as 

f(Ai A„) := F [D 2 u-A{- , u, Du)] , mOT, (2.2) 

where f is a suitably-defined (detailed below), symmetric function of A,, which are the eigenvalues of the 
modified-Hessian matrix: [D 2 u — A(- , u, Du)]. 

In order to start formulating a regularity theory for solutions solving (|2.ip . some conditions need to be 
placed on f (in addition to being a symmetric function), as it is denoted in (|2.ip . First, it is assumed that 
f G C 2 (T) n C°(r) is a symmetric function defined on an open, convex, symmetric region F c R", r/R" 
with F + r + c r, where F + is the positive cone in R". Moreover, it is assumed that f satisfies the following 
conditions: 

f > o in r, f = o on ar, 

f is concave in /", 

and 

5^r}X,-><7! on r^ 2 (0, (2.6) 
/ 

where T Ml M2 = r^ li/l2 (f) := {A e f(f) : /ij < f(A) < /x 2 } for any /x 2 > Mi > and cr , Ci are positive 
constants depending on /xi and /^2. 

Before moving on, some relevant remarks based on the ones made in | SUW04 | are now presented: 
Remarks 2.1. 

(1) Conditions (|2.3p and (|2.4p together imply the degenerate ellipticity condition: 

f/ > 0, in F for / = 1 n. (2.7) 

This combined with the concavity assumption on f subsequently implies that F[D 2 u —A( ■ , u, Du)] is 
a concave function of D 2 u — A{ - , u, Du), which is required to apply the C 2,a estimates of Lieberman 
and Trudinger presented in | LT86 |. 

(2) r(f) enables the definition of an admissible solution to be made corresponding to (|2.ip . That is, a 
solution u G C 2 (fi~) is admissible if 

D 2 u - A{- ,u,Du) e F{f), in IT. (2.8) 

It is clear that (|2.8|) combined with an assumption that 6(x, u) > 0, ensures that (|2.ip is elliptic with 
respect to a solution u G C 2 (Q _ ). In this context, ellipticity and admissibility of functions are indeed 
equivalent and will be used interchangeably for the remainder of this paper. Moreover, a solution to 
(|2.ip is admissible if and only if it is a viscosity solution (see [Tru90j). Indeed, one may substitute the 
term viscosity solution for admissible solution and vice-versa anywhere in this paper. It is through this 
equivalence that the notion of viscosity solutions will be understood for the rest of the forth-coming 
exposition. 



(2.3) 
(2.4) 

(2.5) 
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(3) Assuming (|2.3j) and (\2.4\i to be true, conditions (|2 . 5j) and (|2.6p can be shown to be equivalent to other 
criterion. For instance, if 

f(t t) -> oo, as t -» oo, (2.9) 

then from the concavity of f and the fact that < Yl; it ' s readily calculated that 

f(t < f(\) + ^fi-(t-\i) < f(\) + t^fj. (2.10) 

From this, (|2.5|) follows on any /~o M2 (f) when t is fixed to be large enough. Indeed, (|2.5|) and (|2.9|) 
are equivalent. Moreover, recalling (|2.3p and G dT, (|2.10|) indicates that 

^A,<f(A)<M2, (2.11) 

on r 0:fl2 (f) when taking t = 0. Similarly, if for any and p.2 with < jUi < /x 2 , there is a constant 
9 = 9([Mi, 1X2) such that 

9+f(X)<f(2X), VA G r u,i,^2 • (2.12) 
then the concavity of f can once again be used to directly show 

f(2A) < f(A)+^A,. 

This subsequently implies (|2.5p with a = 9. It is clear that (|2.12|) is satisfied with 9 = (2 a — 1)jUi if 
f is homogeneous of degree a e (0, 1]. 

(4) The main examples of functions f satisfying (|2.3p - (|2.6p are those corresponding to 

f (a) = si /k (\) = ( fn Am )) ( 213 ^ 

\l<ii<—<i k <n \m=ii J J 

and the quotients a^j, 1 <•■•</< k <■■ < n, for which it is denoted that 



S fc (A) 



l/(/c-/) 



As r(Sfc) C T(S/) for I < k, r(f) = r(S^) holds when f is either S^ fc or C/c,/. For these examples, 
the concavity of f is verified in |CNS881 [Tru901 [Urb01 |: and (|2.6p is clear from the argument presented 
in Remark [2~Tl f3")), as both (|2.13p and (|2.14[) depict functions which are both homogeneous of degree 
1. Lastly, in the case where f is defined by either (|2.13p or (|2.14p , it is a straight-forward calculation 
to show that 

f; > 0, in r, for / = 1 n; (2.15) 

i.e. f satisfies the regular ellipticity criterion on T(f) and not just the degenerate ellipticity condition 

(E3 . 

(5) Conditions ()2.3p - ()2.6p combined with (J2.15P are essentially the ones used in jUrbOlj to prove the 
existence of smooth solutions to a class of Hessian equations closely related to (|2.ip . satisfying a 
simpler version of the natural boundary condition depicted in (jl.lbp . 

At this stage, S]/ k and a^,\ represent the two prime candidates for f , for which the existence of globally- 
smooth solutions may be able to be proven. Next, the structure from the Optimal Transportation Problem will 
be used to ascertain criterion on both A( ■ , u, Du) and T u that will subsequently further reduce the possible 
forms of F. 

2.2. Conditions from Optimal Transportation. The next step in the motivation of (|l.la|) - (|l.lb|) , some 
conditions and notation will be recalled from the theory of Optimal Transportation; these will be essential parts 
to the upcoming a priori estimates. Denoting U D Q~ x Q + , the cost-function c (again taking the label from 
Optimal Transportation) is assumed to be real-valued in C 4 (U), satisfying the following three conditions: 

(Al) For any (x, y) e U and (p, q) e D x c(U) x D y c(U), there exists a unique Y = Y(x, p), X = X(q, y), 
such that c x (x, Y) = p, c y (X,y) = q. 
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(A2) For any (x,y) G U, 

Det [D 2 xy c] ^0, 

where D 2 y c is the matrix whose elements at the / th row and j th column is ■ 
(A3w) For any (x, y) e U, and (,r) el" with £ _L 7) such that 

(qy.st - c q,r Cjj, q c rM )c s - k c tJ £jijr) k r)i > 0, 
where Qj(x,y) = 8 q x ^ , and [c' J ] is the inverse matrix of [c,j]. 

Remark. The (A3w) condition is a degenerate form of the following (A3) condition: 

(A3) There exists a constant Co > such that for any (x,y) e U, and £, 7) e R" with £ _L rj such that 

(cy, st - c^c^c^^V^^, > QI6| 2 h| 2 . (2.16) 

The (A3) condition was introduced in |MTW05| as a necessary criterion for local C 2 regularity of the Optimal 
Transportation Equation. However, for the current global regularity analysis, only the weaker (A3w) is needed 
to carry out the necessary estimates provided a technical barrier condition is met (see below). If the barrier 
condition is not met, then the full (A3) condition is needed to ascertain global C 2 regularity. This point will be 
discussed further in Section [B~3l 

The regularity theory presented in |MTW05l ITW06| for the optimal transportation problem gives strong 
indications that A{ ■ , u, Du) and T u need to satisfy certain criterion in order to guarantee that globally smooth 
classical solutions exist in the current scenario. Indeed, much of the structural constructions from optimal trans- 
portation need to be applied to A( ■ , u, Du) and T u if any headway is to be made in proving a global regularity 
result. Specifically, the calculations in Section [5] require that both A( ■ , u, Du) and T u satisfy conditions laid 
down by Optimal Transportation theory in order to make an obliqueness estimate on the boundary condition, 
which plays a fundamental role in the subsequent global C 2 a priori estimate. Moreover, this structure from 
optimal transportation will also be of key importance in the C° estimate presented in Subsection 15.21 

Given the conditions on c, T u is defined as solving 

Du = D x c(- ,T U ), infi", (2.17) 

with A( ■ , u, Du) defined as 

A{- , u, Du) := D 2 c(- ,T U ), in (2.18) 
Thus, the (A3w) condition in the current scenario can be represented as 

D 2 pip A u (x, u, Du)UjVkVi = D pmCu (x,yMjr]kW > 0. (2.19) 

The structure of T u and A( ■ , u, Du) embodied in (|2.17p and (|2.18p respectively, will be assumed along with c 
satisfying (Al), (A2) and (A3w) for the rest of the paper, unless otherwise indicated. 

Remark. Since admissible solutions of (|l.lap - (|l.lb[) must be c-convex, it is readily verified that Loeper's 
counter-example in | Loear | is also valid in the case of (|l.la|) - (|l.lb|) . This — combined with the role played by 
the (A3w) condition in the C 2 estimate reduction to the boundary in Subsection 15.31 — indicates that the (A3w) 
condition is both necessary and sufficient for the existence of globally smooth solutions solving f)l.lafl — (jl.lb[) . 

In the forthcoming obliqueness estimate, the local invertability of the matrix [D 2 u — D x c( ■ , T u )] will be of 
central importance. As this invertability is tantamount to 

Det [D 2 u- D x c(-,T U )] > 0, in fi~, (2.20) 

it is required that T(f) be contained in the cone defined by (|2.20p . Subsequently, the only functions in (|2.13p 
or (|2.14p that have this property are a n j for / < n. In fact, by Remark [27TI| 4|) in Subsection 12.11 r(a n j) is 
equivalent to the cone defined by (|2.20p , when considering (|2.17p and (|2.18p applied to (|2.ip . 

Along with specifying the forms of A( ■ , u Du) and T u , conditions need to be placed on and fi + in order to 
make the forthcoming obliqueness and boundary C 2 estimates. Again, the necessity of these conditions mirrors 
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that of the Optimal Transportation Equation. Specifically, it is required that Q~ and Q + both be bounded and 
C 4 , with ft - , ft + being uniformly c-convex, c*-convex (respectively) relative to each other (See | TW06 | for 
details on c-convexity of domains.). This notion of ft~ being uniformly c-convex relative to ft + is explicitly 
stated as 9ft~ € C 2 with a positive constant 5q such that 

[Dnj(x) - c u c, J , / (x,y)7Ux)]r,T / > 5q , Vx e 3ft", y e ft + , (2.21) 

where r is a unit tangent vector of 9ft~ at x with the outer unit normal 7. If one considers c*(x, y) = c(y, x), 
then (|2.2ip subsequently yields an analogous representation for ft + being c*-convex relative to ft - . This notion 
of ft - , ft + being c-convex, c*-convex (respectively) relative to one another plays a key role in the forthcoming 
obliqueness estimate in Subsection 15. II 

2.3. Conditions on Inhomogeneity. Next, conditions on the inhomogeneity of (|2.1|) will be presented. These 
conditions mirror those stated in | Urb01 | for classes of Hessian equations satisfying a natural boundary condition. 
In order to ensure uniform ellipticity of (|2.ip , it is required that 

e(x,z)>o, xerr.zeM. (2.22) 

In addition, it will also be assumed for all x G ft - that 

£>(x, z) — > 00, as z — > 00, 

B(x,z)-s-0, as z -> -00, (2.23) 

which will be necessary for making the C° solution estimate in Subsection 15.21 Lastly, in order to apply the 
method of continuity, unique solvability of the linearised problem is needed, which subsequently requires that 

6 z (x,z)>0, xeft - ,zeR. (2.24) 

This condition may be relaxed to the following: 

S z (x,z)>0, xer.zet, (|Z24t v) 

via the application of the Leray-Schauder theorem (see | GT011 Theorem 11.6]) as done in | Urb95 |, at the 
expense of uniqueness of an admissible solution. 

Remarks. 

(1) If (|2.22p is relaxed to merely requiring that 6( • , u) be non-negative rather than positive, the eigenvalues 
of the modified-Hessian matrix may not lie in a compact subset of /~(f), even if D 2 u is bounded. Thus, 
it is not possible to deduce the uniform ellipticity of (|2.ip in this scenario. 

(2) (|2.241/vp is also required in the boundary C 2 estimate. 

(3) In general, B(x,z) can not be allowed to have a dependence on Du for the case where f = a n ,i- In 
|Urb95l Section 6], Urbas constructs an example where a solution corresponding to a n j with B having 
a Du dependence has its second derivatives blowing up at the boundary in the two dimensional case; 
this example can be readily adapted to the current class of modified-Hessian equations with f = o n j. 
Indeed, Urbas later applied this example to boundary value problems similar to the ones considered here 
in |Urb97j . Theoretically, in the general case, such a dependence prohibits the existence of barriers 
which are required for boundary gradients estimates subsequently used in the forth-coming obliqueness 
and boundary C 2 estimates. This restriction is again reflected in |Urb01| for a class of Hessian equations 
for the same reason. However, it is possible to have B dependent on Du if the minimum eigenvalue of 
the second fundamental form of 9ft - is large enough for all x e 9ft - . In this scenario, it is possible to 
construct barriers for the subsequent boundary gradient estimates; this is remarked upon in Subsection 
15.11 In addition to this, B also needs to be convex in its gradient argument in order for the C 2 estimate 
in Subsection l5.3l to be applied; this will also be remarked upon further in that subsection. As B having 
any dependence on Du introduces a dependence between B and ft - , it is not a valid criterion for a 
general result. 
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2.4. A Barrier Condition. If A is only assumed to satisfy the weak (A3w) condition, then a technical barrier 
condition is also needed to make the global C 2 a priori estimates in Subsection 15.31 Indeed, if A is such that 
the strong (A3) criterion holds, then no such extraneous condition will be needed; this is discussed in Remark 
15.51 at the end of Subsection 15.31 

Taking A as in (|2.ip , it will be assumed that there exists a function 4> S C 2 (Q~) satisfying 

[D u $(x) - D Pk A u {x, z, p) ■ D k 4>(x)Mj > 5\£\ 2 (2.25) 
for some positive 5 > and for all £ G E", (x, z, p) G U C Q~ x R x K", with pmj n (U) = fi~. 

This condition places a relatively minor restriction on Q~ when /\ is assumed to satisfy the (A3w) condition. 
In | MTW05 |, this condition was needed in order to make the global C 2 estimate for the Optimal Transportation 
Equation; but it was stated in that paper that this condition was removable via a duality argument. For the 
currently considered class of modified-Hessian equations, no such duality exists; thus, (|2.25|) must be kept as a 
separate condition in order to make the following calculations go through. 

2.5. The Quotient Transportation Equation. With the preceding justification regarding the form of F, 
A{ ■ , u, Du) and T u , along with the hypotheses placed on B( ■ , u), attention will now be focused on the following 
set of boundary value problems: 

[D 2 u- D 2 c(-,T U )] = 8(-.u), inn - , (2.25a) 
T U (Q-) = n+, (2.25b) 

where T u is defined by (|2~T7]) : c satisfies conditions (Al), (A2) and (A3w); B(- , u) satisfies (f2~22|) — (f2~24|) : and 
both Q~ and Q + are bounded and C 4 with Q~, Q + c-convex, c*-convex (respectively) relative to each other. 
(|2.25ap - (|2.25b|) will be referred to as the Quotient Transportation Equation for the rest of this paper. 



3. Main Results 

We now are able to state the main results of this paper. 

Theorem 3.1. Let c be a cost-function satisfying hypotheses (Al) and (A2) with two bounded C 4 domains 
Q~ , Q + c M" both uniformly c-convex, c* -convex (respectively) with respect to each other, in addition to either 

• the barrier condition stated in Subsection \2.4\ holding for c and Q~ , with c satisfying the weak (A3w) 
condition or 

• c satisfying the strong (A3) condition. 

IfB is a strictly positive function in C 2 (TF x M) satisfying JZ23]) and (|2.24l/v|) with T u defined by (|2~T7j) . then 
any elliptic solution u e C 3 (Q~) of the second boundary value problem (|2.25ap - (|l.lb[) satisfies the a priori 
estimate 

\D 2 u\<C, (3.1) 

where C depends on c, B, Q~ , Q + and sup^- \u\. 

From the theory of linear elliptic equations, higher regularity automatically follows from better regularity on 
c, Q~, Q + and B. For example, if c, Q + and B are all C°°, then one has that u G C°°(Q~). 

Remark. The dependence of the estimate (|3.ip on sup fi - |u| may be removed if B is independent of u. 

As a consequence of Theorem 13.11 the method of continuity of continuity will be applied in Section [6] to 
prove the existence of classical solutions of (|2.25ap - (]2.25bp . 

Theorem 3.2. If the hypotheses in Theorem \3.1\ hold, then there exists an elliptic solution u G C 3 (Q~) of the 
second boundary value problem (| 1 . 1 a p . (jl.lbp . If, in addition, (|2.24p is satisfied, then the elliptic solution is 
unique. 
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The plan for the proof of Theorem [3.2l is as follows. First, some technical results and inequalities from other 
works will be reviewed in Section |4l From there, various solution estimates will be proven in Section [5l The first 
such estimate will be the obliqueness estimate proven in Subsection 15.11 This calculation will show that the 
boundary condition (|l.lb|) with T u defined by (|2.17p . where the cost-function is assumed to satisfy (Al) and 
(A2), is strictly oblique for functions u where [DT U ] is non-singular. In this obliqueness estimate, the convexity 
assumptions on both Q~ and Q + , play a critical role. This particular estimates differs from the one in | TW06 |, 
as calculations are made without the use of a dual formulation of (|2.25a|) - (|2.25b|) . Instead, the property of T u 
being a local diffeomorphism is exploited to make the argument go through. 

Following the obliqueness estimate, C° bounds for solutions to (jl.lap are then derived in Subsection 15.21 
This calculation mimics the corresponding estimate in | Urb01 | but with a new analogy of a parabolic subsolution 
discovered in | TW06 |. Without this newly-discovered function, it would not be possible to prove the C° estimate 
using the current methods. 

In Subsection l5.3l it is proven that second derivatives of solutions of f|2 . 1|) can be estimated in terms of their 
boundary values, if A( ■ , u, Du) obeys the (A3w) condition ()2.19|) . Again, this estimation deviates from the 
corresponding calculation in |TW06 |, as the case when f = a n j requires the use of specific technical inequalities 
in order to bound the second derivatives in terms of boundary values. This particular argument is carried out 
for general, symmetric A{ • , u, Du) satisfying the (A3w) condition and not necessarily having the specific form 
dictated by (|2.18p . Finally, in Subsection l5.41 the boundary estimate for second derivatives is proven in a similar 
manner to |TW06| : but a key lemma from | Urb01 | is needed to make the estimation go through. This directly 
leads to the global second derivative bounds stated in Theorem 13.11 

Following these a priori estimates, the method of continuity is applied in Section[6]to ascertain the first part 
of the result in Theorem 13.11 regarding the existence of globally smooth solutions uniquely solving (|2.25a|) - 
(|2.25bp . The application of the method of continuity follows the procedure in [TW06]; but instead of simply 
integrating the equation to get solution bounds, the C° estimate from Subsection 15.21 is used instead. From 
there, the Leray-Schauder fixed point theorem is then applied to prove existence of globally smooth (albeit not 
necessarily unique) solutions when one only has B z { ■ , u) > in 

The final section discusses possible directions in research under which to proceed from the current set of 
results. 



4. Preliminary Lemmas 

In this section, some technical relations and inequalities will be recalled from other works that will be subse- 
quently used in the a priori estimates presented in the next section. 



To gain the scope of generality of these lemmas, the following notation will be assumed. Given an arbitrary 
F : R" xn -» R, f : R" ->■ R will be defined by 

f(Ai A„) := F[M], 

where M is an arbitrary n x n matrix having Ai, . . . , X n as its eigenvalues. In addition to this, derivatives on F 
will be denoted by 

F U := °F and F U# d * F 



dMjj ' dMijdMki ' 

with derivatives on f correspondingly written as 

df d 2 f 

Lastly, in the literature, the trace of the operator F is frequently denoted by 

T = T(F) := F" = Y,f- 
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With this notation and that from Section [2] it is now possible to state the supporting lemmas that will be used 
in the forthcoming a priori estimates. 

This first lemma states a technical relation that will be used in the calculations of the C 2 estimate reduction 
to the boundary. This lemma is proven in [Ger96j : but this proof will not be recalled here as it does not aid in 
the understanding of the a priori estimates in Section [5] 

Lemma 4.1. For any n x n symmetric matrix E = [E,y], one has that that 



F ' J,k '-U-kl ~ fl\,.r)X,-"- JJ ' + J2 \[- 1,-U- (4 -^ 



d\id\j'~ JJ ^Xi-X,' 
The second term on the right-hand side is non-positive if f is concave, and is interpreted as a limit ifX; = Xj. 

Remark. Lemma |4~TI was also used in | SUW04 | to make a C 2 a priori estimate similar to the one presented in 
Subsection 15.31 

The following lemma is due to Urbas (see | Urb01 |1: and while the statement of this lemma will be used in 
the boundary C 2 estimate of Subsection [531 the proof itself contains some valuable technical relations that will 
subsequently be used in various locations within the forthcoming a priori estimates. 

Lemma 4.2 ([DrbOl]). If f = a n ,i with I e {1, . . . , n — 1}, then there exists a positive constant C(e) such that 

$>A 2 <(C(e) + e |A|)]r/;, onW 2 (f), (4.2) 



for any e > and < p,i < H2, with C(e) depending only on [L\,\ii and e. 
Proof. Denoting 

S k -ij(X) :=S k -x(X) 

it is calculated that 



= dS k (X) 
a,=o dXi 



n-l \ 5i J V Si Sf 
Next, it is readily observed that 

^S k -j(X) = (n-k)S k (X), (4.4) 

/ 

for k = n. Summing (|4.3|) across / and applying (|4.4p yields 

V- f 1 fSn^ 1 { Sn-^-jn-l+VSnS,-^ 

Next, using the fact that 

5 k (X) = S/ c _i,;(A)A, + S kJ (X) for each / = 1 n, (4.6) 

along with (|4.4|) . it is deduced that 

Sk-MXi = n ■ 5 k (X) - ]T S k .j(X) 

i i 

= k-5 k (X). (4.7) 
On the other hand, (|4.6p can also be used to show that 

5/<-i ; ,(A)A 2 = S k (X)Xi - S k ;i(X)Xi 

i i i 

= S 1 (X)S k (X)-(k + l)S k+1 (X), 
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where 5/< + i(A) is defined to be zero if k = n and (|4.7p has subsequently been used to produce the second 
equality. Using this with (|4.4p and summing, it is next calculated that 

In the special case where / = n — 1, (|4.8p reduces to 



which completes the proof for the case when / = n — 1. 
To proceed further, the Newton inequality is next recalled: 

Sfc(X) S/_i(A) S fc -i(A) 5 ; (A) 

which is valid for any 1 < / < k < n and any A e T + (see | Mit701 Section 2.15]). Taking k = n, Newton's 
inequality yields 

(n-l + l)S n (A)S,_ 1 (A) < ^5 n _i(A)S,(A); 
and therefore, from (|4.5p it is ascertained that 

l/S n (A)\^ S n _ 1 (A) v 1 (S n (\)\^ 5 n _x(A) 

«ls,(A)j S„(A) -2^ r '- „_/ VS,(A)^ S n (A) ■ 1 yj 

Given (|4.8p , (|4.9p and the fact that S„(A)/S/(A) is bounded between two positive constants, the following 
inequality now follows: 

c s, +1 (x) <c Sh-iW , 4im 

for some positive constants Q and C 2 . The ratio S/+i(A)/S n -i(A) is trivially bounded if / = n — 2, so the 
lemma holds in this case as well. 

Now, the remaining cases are considered, which (in light of (|4.10p ) entail bounding the quantity S/+i(A)/S„_i(A) 
from above. Without loss of generality, it is assumed that A = (Ai, . . . , A„) with Ai > • ■ • > A„. Since all A, 
are positive, it is seen that 

k k 

i=i /=i 
for j = 1 n and some positive constants Cj. Therefore, one now has the following: 

5;+i(A) U'tl A/ S;(A) 

- °n5^ Cl ^(A) X/+lAn - C2A,+lV (411) 

Next, it is calculated that 

CoA^ 1 < Ci J] A,- < ^(X) < C 2 ; 

i=i+i 1 

and hence, 

A„ < C. (4.12) 
Now, an arbitrary e > is chosen. If A, +1 < eAi, then by (f4~TTj) and (|4~T2)) 

Jfi-ll.Aj 

If A/ + i > eAi, two cases require consideration. First, if A„ < e, then by (|4.1ip , it is ascertained that 

I^TT ^ CeA '+i ^ CeA i' ( 414 ) 
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while if A n > e, then 

Co*"-'- 1 *! < Ct IT A, < ^(A) < C 2 . 

J/ 



i=/+l 

Thus, one has the following relation: 



Since S„(A)/S;(A) is assumed to be bounded between two positive constants, and S n /S; = on dT + , a positive 
lower bound is now implied: 

C(e) < A n . 

It then follows that 

m^<C (e) . (4.15) 

By applying (|4.13p , (|4.14p and (|4.15p to (|4.10p and replacing e by e/C for a suitably large constant C, it is 
finally derived that 

^•A?<(C( e ) + e|A|)£/) 

for any e > as required. □ 

Remark. The proof of Lemma l4~2l comes straight from | Urb01 |. As mentioned before, this particular exposition 
of proof has been recalled here, as it depicts key calculations and technical relations that will subsequently be 
of use in the forthcoming a priori estimates. 

5. Solution Estimates 

In this section, various a priori estimates to elliptic solutions of (| 1 . 1 a |) . (jl.lbp will be presented. These 
estimates will subsequently be used to prove (via the method of continuity) the existence of globally smooth 
solutions to the Quotient Transportation Equation. 

5.1. Obliqueness Estimate. In this section, it will be proven that the boundary condition (jl.lbp implies a strict 
oblique boundary condition. This estimate will subsequently be used in the continuity estimate in Subsection 
15.21 the boundary C 2 estimate in Subsection 15.41 in addition to justifying the use of the results from | LT86 | 
that yield C 2,a estimates from the forthcoming C 2 a priori bound. 

To begin, a boundary condition of the form 

G(-, u, Du) = 0, on DOT, (5.1) 
for a second order partial differential equation in a domain Or is called oblique if 

G p ■ 7 > 

for all (x, z, p) G dQr x M x W, where 7 denotes the unit outer normal to dQr . 

Next, it is assumed that (p~ and (p + are C 2 defining functions for and Q + respectively; with (p~,(p + < 
near9fi~,<9fi + respectively; (p~ = on dQ~ , + = on 9Q + ; and V$~, V$ + 7^ near9fi~,9fi + respectively. 
If u G C 2 (ft~) is an elliptic solution of the second boundary value problem (|2.25ap — ( |2.25b[) . then the following 
relations hold: 

tP+ o T u = on (p + o T u < near dQ~ . 

By tangential differentiation, it is ascertained that 

4>UDjV)Tj = 0, (5.2) 

for all unit tangent vectors t. Note that the subscript on T u has been dropped without loss of clarity. From 
(|B~2|) . it follows that 

(P+(D,T k ) = xi „ 

for some x > and 7 is again an outer, unit normal of Consequently, one has that 

(t>ic kJ w,, = X7/. (5.3) 
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where 

At this point, it is observed that x > on 9f2~ since |V$ + | ^ on and det DT ^ 0, which is 
subsequently implied by the ellipticity of u and the (A2) condition. Moreover, since u is assumed to be an 
elliptic solution to (|l.la|) , it is observed that 

4>tc kJ = Xw ki lk, (5.4) 
where [w'-i] denotes the inverse matrix of [w-,j]. Upon writing 

G(x,p) := 0+oV(x,p), (5.5) 

one has from the (Al) condition that 

p k :=G Pk (-,Du) = (PlD Pk Y' 

= <t>ic>- k 

= xw k 'li, 

which subsequently indicates that 

P -7 = xw' J irtj 

>0 (5.6) 
on 9Q~. On the other hand, from (|5.4p , it is calculated that 

<t>tc iM Wj k <t>lc iJ = xd>fc iJ 7 J 

= x08-7)- (5 7) 

Eliminating % f rom (|5.6|) and (|5.7|) . yields 

03 • 7) 2 = {w u in j){w kl c^ k c r ^+4>t). (5.8) 

(|5.8p is referred to as a formula of Urbas type, as it was proven in | Urb97 | for the Monge-Ampere equation 
with a natural boundary condition. 

Remark. In the above calculations, the fact that w-,j is invertible has been used; and this is a consequence of u 
being an elliptic solution of (|l.la|) . It is this invertability that requires the equations under consideration to have 
an Optimal Transportation structure combined with a structure that automatically implies that elliptic solutions 
are c-convex (see Appendix A for details). 

Now, P ■ 7 needs to be estimated from below. This calculation mimics the one in [TW06] and |Urb97] for 
the Monge-Ampere equation, but with a modification to avoid the use of a dual formulation to (|l.lap . These 
calculations start with estimating the double normal derivatives of a solution to a Dirichlet problem related to 
(|l.la|) : this follows the key idea from | Tru95 |. Specifically, a point xo G dQ~ is fixed, where p -7 is minimised 
for an elliptic solution u G C 3 (fi~). From there a comparison argument is used to estimate 7 ■ D(p ■ 7) from 
above. Given that P -7 does not have any assumed concavity criterion in the gradient argument, the quantity 
itself needs to be modified so that a workable differential inequality can be derived. Thus, the following auxiliary 
function is defined: 

v := p • 7 — K,((p + o T), 

with a point x on fixed, where p -7 is minimised for an elliptic solution u G C 3 (Q~) for sufficiently large 
k, with the function <p + now chosen so that 

[D u {4> + o T) - c^Qjji ■ , T) ■ D k ((P + o T)] > 5 +|C| 2 (5.9) 

near V£ G R" and some positive constant 5q. Inequality (|5.9p is possible via the uniform c*-convexity of 
Q + with respect to Q~ and taking <p + to be of the form 

cj)+ = a (d + ) 2 - bd + , (5.10) 

where 

d+(y )= mf \y-y \ (5.11) 

yeo£2- 
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with a and b taken to be positive constants, | GT011 |TW06 |. 

Remark. (|5.10p represents only one particular example of </>+ that may be chosen that satisfy (|5.11|) . Indeed, 
there are several alternative barrier functions listed for convex domains in | GT011 Section 14.2]. (|5.10p is a 
generalised construction of a barrier, in that only derivatives up to second order of a barrier work into the 
subsequent boundary gradient estimates. As a bound on the gradient is only needed at it's maximum point 
on the boundary, only a neighbourhood of such a point is considered. In this neighbourhood, a diffeomorphic 
mapping can be applied to straighten the boundary, to see that (|5.10|) indeed represents a second-order Taylor 
expansion about the extremal point x , along the inner-normal of the boundary (see |Eva98 | for examples of 
such calculations). 

For clarity in the forthcoming calculations, the following definition is made: 

H(x, p) := G p (x, p) ■ 7 (x) - kG(x, p), (5.12) 
where G(x, p) is defined by (|5.5P : that is, 

v(x) = H(x, Div(x)). 

Calculating, it is seen that 

D,v = D,H+(D Pk H)u ki 

D u v = DjjH + (D, Pk H)u kj + {D lPk H)u ki + (D Pk H)D k u u 

+ {D PkPl H)u ki u u . (5.13) 

Next, taking F := a n j, equation (jl.laj) in the general form (|2.ip is differentiated to ascertain that 

F ij [D k u u - D k Aij - (D pl Aij)u lk \ = B k + B z D k u. 

Introducing the linearised operator C: 

Cv = F u [D ljV - (D Pk Ajj)D k v)]; (5.14) 
and using (|5.13p and (|5.14p with some simple estimation, it is calculated that 
Cv = F iJ [D u H + 2(D iiPk H)u kJ + {D PkPl H)u kiUlj - (D Pk Aij)D k H] 
+ (D Pk H)(B Xk + B z D k u + F' J D k Ajj) 

< F u D,jH+ F u u, kUjl [D PkPl H + 5 kl ] + F U [(D,, Pk H)(D JiPl H)6 kl - (D Pk A u )D k H\ 

+ (D Pk H)(B Xk + B z D k u + F' J D k Ajj) 

< F' J u lk Uj,[D PkPl H + 5 kl ] + C(F" + 1). (5.15) 

In the above estimation, the gradient bound on the solution u (which is implied by (|2.17p ; the boundedness 
of Q+, and the continuity of c) has been used. Now, with an elementary calculation and condition (Al) 
applied to (|5.5p , one now has the following: 

D PiPj G = D Pj (cP+c k '') 

= (p+c kj c lj -4> + k c 5J c k 5 i 

= c k '> ! c IJ [4>t l -tfc r - s c 5 , kl ] 

Utilising the criterion for (p + in (|5.9p , it is subsequently calculated that 

[D P;Pj G(x,Du)]Uj>5+Y,\ cIJ Zj\ 2 

i 

><5+|£| 2 (5.16) 
for a further positive constant 5^ . Thus, by choosing k sufficiently large, one has 

[D PiPj H{x,Du)]Uj<-\^\ 2 
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holding true near dQ . Substituting this into (|5.15[) yields 

Cv < -^n5 k iF u u ik uji + C(F" + 1) (5.17) 
where C is a constant depending on c, B, ft - , ft + and k. 

Before proceeding, a technicality in the above argument must be addressed. The calculations that lead to 
the inequality in (|5.16|) depend on the explicit structure of cp + . Indeed, unless <p + extends to all of ft + such that 
()5.9|) holds for all T e ft + , there is no control near <9ft~ to validate ()5.16|) and thus (|5.17|) . To get around this, 
one can simply modify the definition of G in (|5.12|) by a function satisfying (|5.16|) in all of ft~ and agreeing 
with (|5.5p near 9ft~. One example of such a G (suggested in | TW06 |1 is given by: 

G(x, p) := p £ * (max{0+ o Y(x, p), C (\p\ 2 - C?)}) , 

where C and C\ are positive constants with Co sufficiently small and Ci > max \Du\ and with e being sufficiently 
small with p € being a standard mollification. 



A suitable barrier is now provided by the uniform c-convexity of Q which implies analogously to the case of 
ft + above, that there exists a defining function <p~ for ft~ satisfying 

[D u d>~ - c''%-,,( • , T)D k (t>-]Uj > 5o\H\ 2 , (5.18) 

near 9ft~. Specifically, <p~ is defined in a similar manner to + in (|5.10|) : 

(p- ■= a(d-) 2 - bd~ , (5.19) 

for some positive constants a and b. Given the elementary fact that D(d~) = —7, the definition above can be 
combined with (|2.21|) and the definition of £ in (|5.14|) to yield 

C4T > F"5o(b - lad') + F iJ _ 2 D,r Dj4T , (5.20) 

in a sufficiently small neighbourhood of dQr; that is, 

Ccp-(xo) > F"5ob+ F'^artrij- (5.21) 

Now the fact that F" is bounded from below on Fq^ 2 is recalled from (|2.5[) and Remark [27Ttf 3]) . With this, 
(|5.17p and (|5.2ip indicate that one can pick a and b in (|5 . 1 9|) such that 

£tfr > Cv, 

in a small enough, fixed neighbourhood of From this, it is inferred by the standard barrier argument (see 

I GT011 Chapter 14]) that 

7 ■ Dv(x ) < C, 

where again C is a constant depending on c, Q + and B. From (|5.18p and since x is a minimum point of 
v on dQr , it can be written that 

Dv(x ) = T 7 (x ) (5.22) 

where r < C. Next, it is calculated that 

D/03.7) = D ; [0+c fcj 7/] 

= cP+D+T 1 )^^ + (Ptc^Dnj + 0+7 J (^ J + c^DiT 1 ) 
= <t> + k c kj {D nj - c s - r ClJ , slr ) + (<j)+ - 0+c r ' s c s , w )c' <J 7 J D,T , 
Multiplying by (p^c tJ and summing over /', one subsequently has the inequality 
<Pfc iJ Dj(l3 ■ 7) = <t> + k <t>+c kJ c tJ {D nj - C^Qj^r) 

+ <Ptc u c k ^ J c l 'dw iq (<P+ - (t>tc r - s c 5:kl ) 
>5 -5>+c''J| 2 , (5.23) 
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by virtue of the uniform c-convexity of Q , the c*-convexity of Q+ and (|5.4p . From f|5 . 12|) and (|5 .22|) , it is 
derived that at xo that the following relation holds: 

-Kw kl c' M c jJ 4>+4>+ < C(j0 ■ 7) - To, 

for positive constants C and To. Hence if j0 -7 < t /2C, one has the lower bound 

w kl c'' k c^<P+(t>+ > ^. (5.24) 

At this point in the argument, a Legendre transform (or in the case of Optimal Transportation a c-transform) 
would usually be invoked to derive a dual equation to (jl.lap . with the intention to bound w/^7,7/ at xo from below. 
In the current situation, such a dual formalism can not be explicitly constructed; indeed, the transformed Quo- 
tient Transportation Equation would not have the same form as (|2.25ap - (|2.25b[) . Instead of doing this, (|5.22p 
is manipulated another way via the fact that T u is a local diffeomorphism (In fact, it's a global diffeomorphism.). 

At x , ([B~22]) indicates that 

w' Jr yiDjV = TH/ y 7/7j 

< CV 7 7,7/. (5.25) 
Using (|5.3p and following the same vein as the calculation done in (|5.23p , one calculates that 
w'hDjiP • 7) = x-^tc^Diitfc'^j) 

= x - 1 <t>t<i>ic k ''c , ' J {Dn J -c s - t c u , 5lt ) 

+ X - 1 4>tc kJ w iq c l -<i 1 jc kJ (<l)+ - 4>U ns Cs, k i) 



><5 +El7/<c u | 2 , (5.26) 
/ 

where 5^ and 5q are the constants associated with the uniform c-convexity of Q~ and the c*-convexity of Q + 
respectively. 

Combining (|5.25P with (|5.26p and using the definition for v along with (|5.7|) . it is ascertained that 

w u r,D jV > 5+ E \lkC kJ \ 2 - x-^4>i^C' k c jJ w kl 
1 

= <5 +ElT/<c w | 2 -4/3-7), 
at x . Setting t + := 5q ^2 k , \'jkC k ' l \ 2 at xo, the above reduces to 

w u r,lj> ^r-^(/3-7). 



Thus, if (/3 ■ 7) < then it follows that 



w U inj > (5.27) 

Remark. In order to estimate (pf^c'^c^'wki, only the positivity of the second term of the right-hand side of 
(|5.23p (due to c*-convexity) was used. That is, the estimate for (pf^c'^c^'w^i only depends on the value of 
5q and not-necessarily uniform c*-convexity of Q + . On the other hand, the estimate for w'-'jijj depends on 
5q and the not-necessarily uniform c-convexity of ft - . 
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Combining (|5.8p , (|5.24p and (|5.27p , the uniform obliqueness estimate is thus derived: 

G p -7><5, (5.28) 

for some positive constant 5 which depends on c and B. This estimate is now stated formally in the 

following result. 

Theorem 5.1. Let c £ C 3 (IR x R") be a cost-function satisfying conditions (Al) and (A2) with respect to 
bounded C 3 domains Q~, Q + £ K", which are respectively, uniformly c-convex, c* -convex with respect to each 
other. In this case, if B(x,z) is a positive, bounded function in C 1 (Q _ x M") and T u a mapping satisfying 
f|2. 17|) . then any elliptic solution u £ C 3 (Q~) of the second boundary value problem (|2.25ap -( |2.25b|) satisfies 
the obliqueness estimate (|5.28p . 

Remarks 5.2. 

(1) The key difference in this obliqueness proof is that a dual formulation to (|l.la|) is not used anywhere 
in the estimate. Instead, the fact that T u is a local diffeomorphism is exploited to make the estimate 
(fS~24)) and (f5~27|) go through. 

(2) As mentioned earlier, the restriction that B be independent of the gradient in the general case is 
unattractive. However, a similar criterion was required by Urbas in | Urb01 | in the case where equations 
involved the operator o n j applied to Hessian matrices. The reason for this restriction is due to (|5.20p 
holding for general barriers and the fact that F" can remain bounded even as \D 2 u\ — > oo. This 
observation is clear, recalling the relation (|4.5p derived in Section|4]and applying MacLauren's Inequality. 

(3) This estimate can be derived in special cases where B has a Du dependence provided that the second 
fundamental form of dQ~ has a large enough minimum eigenvalue. Specifically, it is required that 

2/ r ' > [G,, PkPs J k - kG Kp ,\ + {G PkPrlk - KG Pr )B Ps < -2F ir G PkPs D nk . (5.29) 

This can be realized by examining (|5.15p and noting the appearance of a H Pr B Ps u rs term when B has a 
Du dependence. If (|5.29p holds, this new term can be removed immediately from the calculation; and 
thus, it would not affect the subsequent barrier construction. As f|5.29|) implies a dependence between 
B and Q~, this situation can not be assumed for the general case. 

5.2. C° Estimate. In [MTW05], C° solution estimates are naturally handled as the Optimal Transportation 
problem corresponds to a Monge-Ampere equation, which subsequently correlates to the determinant of a 
Jacobian corresponding to a measure-preserving coordinate transform. Thus, simply integrating the equation 
will lead directly to uniform integral bounds that, in turn, directly lead to C° estimates of the solution. In the 
case of the Quotient Transportation Equation, there are no such simplification given that <r n j has no direct 
correspondence to a Jacobian of a non-singular coordinate transform. Instead, the estimate set forth by Urbas 
in | Urb01 | will be adapted to the current situation, with the use of a special auxiliary function first introduced 
in ITW06I . 

5.2.1. Supremum Solution Bound. First, an upper-bound for an elliptic solution u of (|2.25ap -( |2.25bp will be 
proven. The argument for this part of the estimate is the same as the one |Urb01| and needs a trivial modification 
for the current scenario. Recalling (|2.10p . the concavity of f implies 

fW<f(l + 

and thus, 

F[u] < Ci + aA x (u- c), (5.30) 
where a = JT fj(l, . . . , 1). Integrating (|5.30p and applying the divergence theorem, it is readily calculated that 

/ B(x, u) dx < Ci Vol(ft-) + a [ A x (u - c) dx 
Jq- Jq- 

< C 2 Vol(Q-) +cr / (j-Du)dx 

Jdn- 

< C 3 , (5.31) 
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where 7 is the unit outer-normal vectorfield relative to dQ~ . The last inequality comes from the gradient 
solution bound implied by (|2.17p : the continuity of c; and the boundedness of both ft~ and ft + . From (|5.3ip 
and (|2.23p , it is clear that u is bounded from above somewhere in ft~; thus, by the aforementioned bound on 
Du, it is ascertained that 

sup u < C, 

Or 

where C depends on B, c, ft - , and ft + . 

5.2.2. Infimum Solution Bound. Before proceeding to the infimum estimate, the following definition is needed. 

Definition 5.3. A c-convex function (p on an arbitrary domain U is said to be uniformly c-convex on U, if (p 
satisfies the following inequality 

[<Pu ~ Cij( ■ . T u )} > 0, in U, (5.32) 

where T u is defined by (|2.17p . 

Remark. It is clear that in the case where f := c n j, (|5.32p corresponds to the ellipticity criterion (|2.8p . 

To get the lower bound for u, a specific type of auxiliary function ivo(x) is considered, which is uniformly 
c-convex. This auxiliary function was first used as the basis for an alternate proof of Lemma 5.1 in | TW06 |. 
Let y be a point in ft + and u be the c*-transform of the function 



ip(y) = -Vr 2 - W - yo| 2 , 

given by 

u (x) = sup [c(x,y) -ip(y)], (5.33) 

yes r (y ) 

for sufficiently small r > 0. u is a locally uniformly c-convex function defined in some ball B R (0), with R — > 00 
as r — > 0; and the image of its c-normal mapping satisfies 

T Uo (Q-) c S r (y ), 

where T Uo is a diffeomorphism between Q~ and T Uo (Q~). The reader is referred to | TW06I ITWar | for the 
derivation and discussion of these properties of u . 

With the above construction, it is clear that for given x e Q~ and yo € Q + , r > can be fixed small enough 
so that one has u is uniformly c-convex on fi~ <s B R (0) and 

T Uo (Q-) c e r (y ) CQ+. (5.34) 

Since, Br(0) <s Q _ (by the smoothness of both u and c combined with (|5.32p ), it is clear that 

'V 
Si, 

for some positive constant C(r). Next, it is supposed that 

B{-,u)<C{r), \nQr. (5.35) 

From this, one has that 

j) " ' [D 2 u - D 2 x c( ■ , T U0 )} >B(-,u)= "' [D 2 u - D 2 x c( ■ , T u )], 

so uq — u obtains its maximum on this maximal point is denoted xo. With that, it is calculated that 

D 7 -(u - u)(x ) > and 5j(u - u)(x ) = 0, 

where 7~ is the unit normal vector of Q~ at xo and denotes the tangential part of the gradient; that is, 
5j := D — 7(7" • D). By a suitable rotation of coordinate axes, it is assumed that 7~ = ei at xo; so one now 
has the following: 

DiUo(xo) > Diu(xo), 

D a u (x ) = D a u(x ), for a = 2, . . . , n. (5.36) 



[D 2 u - D 2 c(-,T LI0 )} > C(r), in ft" 
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Using the definition of T u , explicitly stated in (|2.17p , the above is seen to be equivalent to 

ci(x , T Uo (x )) > ci(x , Tu(x )), 

c a (x , T Uo (x )) = c a (x , T u (xq)), for a = 2 n. (5.37) 

From here, it is recalled that [DT U ] > 0, as u is an elliptic solution to (|2.25ap - ()2.25b|) with c satisfying condition 
(A2). Using this fact and the obliqueness estimate (|5.28p , it is seen that the pull-back of 7 + , the unit normal 
vector Q + at T u (x ), satisfies 

< i+DiT k nr 

= 1 +w IJ d' k 1 r. (5.38) 

Since Wjj is positive definite, as u is an elliptic solution of (|2.25a|) — (|2.25b|) with [DT U ] > 0, (|5.38|) along with 
d- k being diagonalizable, as it is symmetric in its indices, implies that 

ci.i(*b.T u (xb))>0, (5.39) 

with our coordinates again labeled such that 7~ = ei. Moreover, (|5.38p implies that po = T u (xq) must lie in the 
set dQ + := {p G dQ + : 7^~(p) > 0}. Thus, the y-gradient of c x at (xo,Tu(xo)) has a non-trivial component 
in the normal direction of at xo. This combined with (|5.37p subsequently implies that 

("r uo (xo)-T u (x ))-7" >0. (5.40) 

However, by taking c*(x,y) = c(x,y) in (|2.2ip and pulling back by the mappings c x (x, •), it is seen that 
being c*-convex relative to fi + implies that the image of c x (x, fi + ) is uniformly convex with respect to x e 
|TW06 |. Subsequently, combining this observation with (|5.37p and (|5.40p , thus implies that 7" Uo (xo) lies outside 
of Q + , which contradicts (|5.34p . Therefore, (|5.35p is false, and so it must be that B(x, u) > u(r) somewhere 
in Q~. Thus, the bound on the gradient of u implies a lower bound on the solution: 

C < inf u. 

Or 

With this, the following result is thus proven: 

Theorem 5.4. Let c e C 3 (R x R") be a cost-function satisfying conditions (Al) and (A2) with respect to 
bounded C 3 domains Q~ , Q + € R", which are respectively, uniformly c-convex, c* -convex with respect to each 
other. In this case, ifB(x, z) is a positive function in C 1 (Q _ x R") satisfying (|2.23p and T u a mapping satisfying 
(|2.17p . then any elliptic solution u e C 3 (Q _ ) of the second boundary value problem ()2.25ap -( )2.25bp satisfies 
the solution bound: 

sup|u|<C, (5.41) 

n- 

where C depends on c, B, fi~ and Q + . 
Remarks. 

(1) In the infimum estimate the obliqueness estimate was used, which subsequently utilised the fact that 
Y^,i fi is bounded away from zero on T j i(i2 . It is important to note that the obliqueness estimate does 
not depend on n\. Obviously, if this were the case then the infimum estimate above would fall victim 
to circular logical argument. 

(2) Applying (|5.4ip to the condition (|2.23p , indicates that B is bounded from above and away from zero. 
Thus, the forthcoming calculations pertaining to (|2.25a|) — (|2.25b|) will be carried out on /"^^((Tn,/), 
where < p,\ < p.2 without loss of generality. 

5.3. Global Second Derivative Estimates. In this section, global bounds for second derivatives of elliptic 
solutions of equation (| 1 . la[) will be shown to be estimated in terms of their boundary values. The following ar- 
guments mimic those presented in [TW06] and [MTW05| . which are subsequent modification of the arguments 
presented in [GT01I Section 17.6]. 
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For this estimate, it suffices to consider the more general form of the Quotient Transportation Equation 
depicted in (|2.ip , under the assumption that the matrix valued function A £ C 4 (Q~ x K x R") satisfies 
condition (A3w); that is 

D 2 pkPi Aij(x,z,pMi7] k Vi > 0, 
for all (x,z,p) £ fi~ x K x K", £,?7 £ R" and £ _L 77. It is also assumed that is /A is symmetric, hence 
diagonalisable. Lastly, the barrier condition from Subsection 12.41 is also assumed to hold for A. This condition 
is recalled here for convenience: it is assumed that there exists a function <p £ C 2 (ft~) satisfying 

[D u $(x) - D Pk A u {x, z, p) ■ D k $(x)]Uj > ~5\tf. (5-42) 

for some positive 5 > and for all £ £ K", (x,z,p) elic^xMxI", with proj n -(i7) = Remark [531 
at the end of this section, will discuss how this barrier condition can be removed if A satisfies the strong (A3) 
criterion. 

To begin, u £ C 4 (fi~) is assumed to be an elliptic solution of (|2.ip with (x, u(x), Du(x)) £ U for x £ 
and £ £ S". The auxiliary function v is now defined as 

v := \oq{Wijiiij) + K(p, 

where w t j := u,j — A-,j. Differentiation (|2.ip yields 

F ij [D^u u - D^A U - {D z A u )D iU - (D Pl Aij)u^] = B £ + B z D^u. 

Another differention subsequently produces 

F^'D^jD^Wk, + F ij [DtfUij - DftAij - {D zz A u ){D^u) 2 - {D PkPl A u )u ik uy 

-(D Z A>« - {D Pl A u )D iUli - 2{D tz A u )D iU - 2(D iPk A u )u ik 

-2{D zPk A u )u ik D iU \ = D ££ B + (D zz B)(D iU ) 2 + 2{D iz B)D t u + {D z B)u ti . 

Recalling that a gradient bound for an elliptic solution solution u exists (implied by the boundary condition 
(jl.lbp and (|2.17p ), one can write the above as 

F'i ,kl D^WjjD^wu + F' J [D&u u - (D PkPl A u )u ik u v - {D Pk A u )D k u i( \ 

>-C[l + (l + WjJ )F"] 

> _ C [(1 + WjJ )F"] , (5.43) 

where the last inequality comes from the fact that < a < F" on Fo iM2 . To further reduce this, the (A3w) 
condition is now used. Specifically, a point xo £ Q~ is fixed and coordinates are rotated such that [wjj] (and 
thus F' J ) is diagonal. With such a rotation, it can also be enforced that the eigenvalues of [w,j\ are ordered as 
< A n < ■ • • < Ai. Estimating, one sees that 

F iJ D PkPl AijD ki uDi C u = F ij A ilkl {w ki + A ki ){w ti + A ti ) 

> F u A Uik ,w^w, i - C(l + Wjj)F n 

> frArr, k l(^ k )(Ml)-C(l + Wjj)F" 

k or l—r 

> -C(F" + w JJ F ii + \w\fj\i) 

> -C(F H + \w\ + \w\F"), (5.44) 

where \ w\ = J2i w a'< ar, d the f act that fiA, is bounded on F MliM2 (as shown in f)2. 1 ip ~) has been used to ascertain 
the last inequality. Using this, (|5.43p can now be rewritten as 

F«' [D ljUii - (D Pk A u )D k u^] > -F* k 'Dew u Dew k , - C(F" + w„ + w n F"). (5.45) 

(|5.45p is used to guide the definition of the linear operator that will be used to analyse v. 

Cux := F'J [D ljUei - (D Pk Au)D k u^ . 
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From (|5.44p , it is calculated that (recalling w,j = u,j — Ajj) 

Cwtf > -F U ' k 'D^w u D i w k i - C(F" + wu + w u F u ), 
for a further constant C. Next, the auxiliary function is differentiated to get 

DjV = - ! —^- + KDj<t> 

DijV = Dum _ ^mDjm + KDiil 



w, 



(5.46) 



ft 



Given (|5.42p , it is clear that C<j> > F"5. Using this fact along with (|5.46p , the following inequality emerges: 



1 



ft 

>- — F^D iWl] D iWkl -F<> DiW ^ W tt 



w. 



ft 



— (F" + wii + WnF") 



+ K.5F". 



(5.47) 



Next, it is supposed that v takes its maximum point at x € in a direction corresponding to the vector £. 
Coordinates are then relabelled so that £ = e±. To proceed, the first two terms on the final line of (|5.47p need 
to be estimated at this maximum point. Specifically, it is claimed that 



-^F^D lWlJ D lWkl -F^ D ' W ^^ >0 
wu wU 



that is, 



— wuF IJ ' DxWjjDiWki > F' J DjWuDjWix. 
To show (|5.48p . it is first claimed that the following inequality holds: 



Wll p>jMi DiWjtDiWki > — Xi 



|D r w n | 2 



r=2 



(5.48) 
(5.49) 

(5.50) 



This relation will be fully justified at the end of the proof. With (|5.50p . one is able to bound the left-hand side 
of (|5.49p from below. Using this, it is observed that (|5.49p is true if 



Ai 



Ai — A r 



> 



J2fr\Dr 



wn\ 



holds, where coordinates have been rotated so that F' J is diagonal. By breaking this inequality up term-by-term 
and dividing each resulting inequality by |D r wn| 2 , one sees that if 



-2Ai - X ~. r > f r , for re {2 n} 

Ai — A r 



(5.51) 
(5.52) 



are shown to be true, then (|5.50p and hence (|5.47p will be proven. For the following calculations, coordinates 
are relabelled so that Ai > A2 > • ■ • > A„, keeping in mind that A„ > jUi > for A e /~ MlM2 . (|5.5ip will first be 
verified. To begin, (|4.3p is differentiated to ascertain 



d 2 f 
9AT 



n- I J \ n - I 

2 



^-2 



n- I 



S;S n -l;, — S/_1;/S n 

s 2 

S n \ ~ 1 ( S 2 _ 1 . ; S n — S|S/_i;,'S n -i;/ 
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Using this with (|4.3|) along with (|4.4|) . it is calculated that 



Ai d 2 f 2AiS/_i,i /" 1 \ / SjS n -i ; i — S n S/_i ; i . ^ 



fi S, \n-l ) \ S n S 

2AiS/_i ; i ( 1 i i ' 



5/ 



Si \n — I J \ Si 

2AlS/-l;l + S/;l 1 S/-1 



l + ^[XiS,- 1;1 -^-) . (5.53) 



S; n - I Si 

. AlS/_l;l ~ 

S, V n- I, 

Next, the observation is made that 

S, ; i < (n-/)A 2 S,_ 1;1 , /e{l n-1}, 

with equality holding when all A, equal and I = n — 1. Thus, it is seen from (|5.53p that 

"XaAf - 1 + J, [ XlSl - 1 ' 1 —/ ) 

= 1 + %^(A!-A 2 ) 

> i, 

which directly implies (|5.5ip as f, > on /~ MliM2 . 

Next, (|5.52p is proven for arbitrary re{2 n}. Again, using (|4.3p and (|4.4p , one finds the following: 

Aifi _ Ai(5 n -i;i5; — 5 n 5;-i ; i) 

Aj/J \j(S n -ijSi — SnSl-ij) 

(S n — S n ;l)Sl — (S; — S/;i)S n 

(S„ — S n j)Si — (Si — Sij)S n 

Sl-lSn — Sn-iSi 
SljSn SnjSl 

Si^i 
Si-j 

AjS/_l;lJ + S/;1J 



AlS/-l;lJ + S/;1J 



(5.54) 



where the fact that S n; ,- = for any / has been used to gain the fourth equality. From the current selection of 
coordinates, it is clear that the following inequalities hold: 

2X] <Xl + AiAj, 
2Xj < Xi + Xj. 



This combined with (|5.54p yields 



2A, ^ Ai/j < 2A J 2 5;_i : i J + 2A i 5/ : i, J 



Xi+Xj Xjfj - (A? + XxX^Ss-x-u + (Ai + Xj)Si; ld 
< 1, 

which subsequently implies (J5.52P via an elementary algebraic manipulation, using the fact that f, > on /~ M 
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In light of the above calculation, the claim in (|5.48|) has been proven. With this, (|5.47[) is further reduced 
to find that 

Cv > — — ( F" + w,j + WjjF'") + kF" 



> -C 



F" 

1 + + F" 

wu 



kF' 



> -C— + (k - C)F", (5.55) 

Wn 

where the fact from Remark [2~TT |2]) that F" > a > on /~o iM , for any /z 2 has been used to gain the last 
inequality. On the other hand, since x is a maximum of v, one has that 

Cv(x ) = F u [D uV (x ) - (D Pk Aij)D k v(xo)] 

= F iJ D ijV (xo) 

< 0. 

This, combined with (|5.55p yields 

> w u (xq), (5.56) 



K — C 

provided k > C, where C depends on sup \u\, sup \Du\, A, sup \B\, Q~ and Q + . 

With this, all that is left to be done to prove (|5.56p is to demonstrate that (|5.50p holds in the case where 
f = (T n j. To do this, it is first calculated that when f = a n j, one has that 

d 2 f ( 1 ,"\ ( 1 \ ( S n \"-' 2 ( S n -ij S n Si-ij\ ( S n -i ; / S n S;_i ; , 



1 



u dXjdXi \n-l J\n-lJ\S,J \ S, Sf J \ S, Sf 

— — l 

1 (Sri]"' 1 f S n -2;ij S n -i-iSi-i-j Sn-l-jSi-i-j S n S/_2;/j ^ ^nSl-l-jSj-i-j 



n-l \ S,J V 5, Sf Sf Sf Sf 

(5.57) 

Next, applying the ellipticity criterion of f, > from (|2.15p to (|4.3p one sees that 

S n Si-ij < S n _i;,S;. (5.58) 

From (|5.58p it is clear that the first term on the right-hand side of (|5.57p is negative. To analyse the second 
term on the right-hand side of (|5.57p , it is first calculated that 

S n -2;ij S„Si-2;ij _ S n S n S/_2;/j 

Si Sf ~ XiXjSi Sf 

S n A XjXjSl-2;ij 



XiXjSi \ Si 
S n ( ^ Si — Si-jj 



XiXjSi \ Si 

SnSlJj 



(5.59) 



XiXjSf 

This enables one to write the quantity inside the parentheses contained in the second term on the right-hand 
side of (|5.57p as 

S n Sl-ij (S n — S n ;i)(Si — Slj) (S n ~ S n -j)(S/ — S/ ; j) ^S n (Si — S/;/)(S/ — Slj) 

XjXj XjXjSf XjXjSf sf 

= %F ( s ' s ';'7 " s ' s ';< " s ' s ';j + 25/ ; /S, ;i ) (5.60) 
■->/ 

Using the observation that 

Si = Sij + Si-j - Si-jj, (5.61) 
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on (|5.60p enables one to conclude that quantity inside the parentheses contained in the second term on the 
right-hand side of (|5.57p is equal to 



S n - 



U [S W (S U + S,j) - (S, 2 , + Sf u )] < 0, (5.62) 



s 3 



where the inequality is immediately obvious. Thus, it has been shown that 

8 2 f 



dXidXj 

i.e. o n i is indeed concave. 



<0; (5.63) 



Defining M,j = D\WnD\Wjj where w corresponds to the w in Section IB~3l and rotating coordinates so that 
Mjj is diagonal with wu corresponding to the maximum of D 2 w at xo, it is clear that the second derivative 
terms on the right-hand side of (|5.50p are all individually non-negative. 

Next, the quotient terms in (|5.50p are analysed. 

fi — fj 1 1 I S n \ "-' ( Sn-ljSl — Sn-l\jSl — SnSl-iJ + S n S/-lJ 



Xi - Ay Xj - Xj n - I \ Si J \ Sf 

1 1 fSn^- 1 S n 



x, 




A, 


n-l 




1 




1 


Xi 




A/ 


n-l 




1 




1 



Si J X,XjSf 



[S, (Xj - Xi) + XiXj (Si-u - S;_i ; /)] 
;2 [S, (Xj - Xi) + X,(S, - S,.j) - Xj(S, - S,;i)} 

Xi-Xjn-l\W^^ [X ^ ! - XiSljl (564) 
It is clear that A, — Ay will have the opposite sign of XjS/j — A,S; ;J . Thus, it is seen that the quotient terms 
in (|5.50p are individually positive, which is a stronger statement compared to the case when f is just concave. 
Indeed, it is the specific structure of cr n j that enables one to discard the terms that result in (|5.50p . With this, 
(|5.50p has been shown to hold in the case where f = a n j, hence proving (|5.56p . 

It is recalled that (|5.56p has been shown under the assumption that v obtains it's maximum at an interior 
point of Q~ . If v does not take a maximum at an interior point, then it must take a maximum on dQr . With 
this observation, the reduction of the C 2 estimate to the boundary is now complete. 

Remark 5.5. In the case where A satisfies the strong (A3) condition, the above estimate is much simpler; 
and does not require the barrier condition (|5.42p . Specifically, one only needs to consider the simpler auxiliary 
function v := in this case, with the definition of £ remaining unaltered. Upon noting that the (A3) condition 
directly implies that 

F'JD PkPl A u D ki uD li u>CF"\D 2 u\ 2 , 
for some positive constant C, it is a straight- forward calculation that shows 

Cv > F'"(Cx\D 2 u\ 2 - C 2 (\D 2 u\ + 1)). 
From this, the second derivative estimate immediately falls out as Cv < 0. 
With this remark, the following has now been proven: 

Theorem 5.6. Let u be an elliptic solution to (|2.ip in Q~ , with x, u(x), Du(x) G U for all x G Q~ and U 
bounded. If (|5.42p is true, with B(x, u) a positive, bounded function in C 2 (Q~ x K), with either 

• A satisfying the weak (A3w) condition plus A andQ~ satisfying the barrier condition stated in Subsection 
Elor 

• A satisfying the strong (A3) condition; 



SECOND BVP FOR MODI Fl ED-HESSIANS EQUATIONS 



23 



then one has the following estimate: 



sup|D 2 u| < C 1 +sup|D 2 u| ) , (5.65) 



where C depends on sup \u\, sup \Du\, A, sup £>|, Q and fi + . 

From these calculations, it is seen that the (A3w) criterion for A on the set U is indeed sufficient to make 
the argument go through. Also, it is clear that the calculation holds for the Quotient Transportation Equation 
(|l.la|) : thus, the calculations of this subsection provide a crucial step to the proof of Theorem |3. II 



Remark. It is possible to generalise Theorem 15. 6l to the case where B(x,z,p) > with 

B PkPl (x.z,p)tkti>0, (5.66) 

for all (x, z, p) £ U. To do this generalisation, (|5.42p needs to be strengthened. Specifically, it is required that 
A and B be such that there exists a 4> satisfying 

[D u $(x) - D Pk (A,j(x, z, p) + B(x, z, p)) • D k $(x)]Uj > 5|£| 2 , 

for all (x, z, p) e U. With that, the calculations proceed precisely as above with the linear operator C now being 
defined as 

Cu := F ,J [D u u - (D Pk A,j)D k u} - (D Pk B)D k u. 
This combined with Remark l5~2Tf 3|). indicates that the only way for B to have a dependence on the gradient of 
u and still allow for the existence of globally smooth solutions to the Quotient Transportation Equation, is if B 
satisfies (|5.66p with Q~ of sufficiently high normal curvature so that (|5.29p holds. 

5.4. Boundary estimates for second derivatives. In this subsection, the C 2 a priori bound is completed by 
proving second derivative bounds for solutions of the Quotient Transportation Equation on dQ~ . This treatment 
is similar to the one presented in | MTW05llLTU86IIUrb98 |, but requires some modification to accommodate the 
particular situation where f = a n j. Specifically, Lemma [4~2l is required in order to derive differential inequalities 
from which barriers can be constructed. A key point to the following argument is that for oblique boundary 
conditions of the form (|5.ip , where the function G is uniformly convex in the gradient, the twice tangential dif- 
ferentiation of (|5.ip yields quadratic terms in second derivatives which compensate for the deviation of /3 = G p 
from the geometric normal. 

First, non-tangential second derivatives are treated. Taking ^ e C 2 (Q~ x R), the following definition is 
made: 

v := M/(x, Du). 

Defining the linear operator C by (|5.14p and calculating as was done for (|5.15p , one observes that 

\Cv\ < C 1 F iJ 5 kl u ik Uj l + C 2 (F ii + 1), 

< C M + J2 fi&f + !)1 ■ ( 5 - 67 ) 

In the last inequality coordinates have been rotated so that F' J is diagonal, without loss of generality. Using 
Lemma l4~2l and defining 

M := sup |D 2 (u - c)|, 

q- 

(|5.67p is further reduced to 

\Cv\ < (C(e) + eMF") + C(l + F"), 

< C(C(e) + eMF") 

< (C(e) + eM)F", 

where the fact that F" is bounded away from zero has been used along with a rescaling of e by a factor of C, 
again without a loss of generality. 
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Using the construction of <p + from Subsection 15.41 is now set in the following manner: 

U/(x, Du) = G(x, Du) = (p + o T u (x, Du). (5.68) 

(|5.19p can be used to construct both an upper and lower barriers (using different choices of a and b), as (|5.68p 
indicates v = on dQ~ . By using the same barrier argument from the obliqueness estimate, the following 
boundary estimate is derived: 

\DG\ < (C(e) + eM), on 
for any e > 0. This, in turn, implies that 

D(/3- Du)\ < C(e) + eM, on dQ~ , 
where it is recalled that P := G Pk . With this, upp is thus estimated on the boundary, which is equivalent to 

Wpp < C(e) + eM, on dQ~ . (5.69) 
Remark. By the strict obliqueness estimate (|5.28|) , it is clear that 

w 77 < Cwpp, (5.70) 

where 7 is the outer unit normal to dQ~. 

To proceed, an explicit representation of any given vector £ is written in terms of a tangential component, 
t(0, and j0: 

e = T(£) + |^A (5.71) 

where 

r(0 = $-(^-7h-|^/3 T , 

and 

/3 T = j0 - 09 • 7)7. 

From this, it is calculated that 

weP = i-4£e--r)-(i-|^ 

<l-2^K-7) + CK-7) 2 . 

It is now assumed that the maximal tangential second derivative of w over dQ.~ occurs at a point at x e dQ~ 
in the direction which is taken to be ei. Denoting r = r(ei) and utilising (|5.7ip , it is calculated that 



w n = w TT + — z w tP + (pt 1 wpp. (5.72) 



27i 7i 2 
Next, the boundary condition 

G(x,Du) = 0, ondn~, 
is differentiated in the tangential direction to find that 

= G Xt +TjG Pk u ik 
= G Xt + u Pt ; (5.73) 

that is, 

Wp T < -G Xr + C 

< c. 
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Subsequently, this combined with (|5.72[) and (|5.69p yields 

w u < |r| 2 Wll (x ) + cf^ + (C(e) + e/W) 7l 2 



< (1 - 2^(£ • 7) + C« ■ 7)^J HtaW + 

+ (C(e) + eM) 7 2 . (5.74) 

Now the term DjgWn(xo) will be estimated. To do this, (|2.1|) is differentiated twice to find that 

F u D u Wu = -F'^'DxWjjDxWk, + B u + 2(B z u u + B lz Dm) + B zz (D 1 u) 2 
> C, 

where the concavity of F has been used along with the assumed condition of B z > 0. Using this and the barrier 
construction in |GT01l Corollary 14.5], with a linear operator defined by 

Cu := F u DjjU, 

one obtains 

- Dp Wll (xo) < C. (5.75) 
Next, (|5.73p is differentiated again in the tangential direction to get 

Gp k p,UkTUiT + Gp k UkTT + G XtXt + 2G XrPk Uk T = 0. 
This combined with (|5.75|) and {|5 . 16|) indicates that at xo one has that 

CoiVn(xo) - CiUn(xo) < C 3 , 
for positive constants Co, C\ and C z . That is, 

wu(xo) < C 

Using this relation in (|5.74p thus yields 

w n <C(e) + eM, on dQ~ , 
for new constant C(e) and a rescaling of e. 

Combining this last relation with (|5.69p and (|5.70p allows one to bound the second derivatives of u on the 
boundary of Q~ in the following manner: 

sup |D 2 u| < C(e) + eM. 

8Q- 

Now utilising (|5.65p to eliminate M from the above inequality finally yields 

sup|D 2 u|<-^, 
an- l - ec 

for a new C, C(e) and e scaling. Fixing < e < the following has been proven: 

Theorem 5.7. Let c be a cost-function satisfying hypotheses (Al), (A2), (A3w) with respect to bounded C 4 
domains Q.~ , fi + £ R" which are respectively uniformly c-convex,c* -convex with respect to each other. Let B 
be a strictly positive function in C 2 (Q~ x M x R") satisfying (j2.23[) and (|2.24l/v|) . Then any elliptic solution 
u e C 3 (Q~) of the second boundary value problem (jl.lap . (jl.lbp satisfies the a priori estimate 



sup|D 2 u|<C, (5.76) 

an- 



wnere C depends on c, B,Q , fi + and sup^- u 
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This combined with Theorem 15.61 proves the main result Theorem 13.11 Once the second derivatives are 
bound, (| 1 . 1 a p is effectively uniformly elliptic. This combined with the obliqueness estimate yields global C 2,a 
estimates for solution of (|2.25ap - (|2.25b[) , from the theory of oblique boundary value problems for uniformly 
elliptic equations presented in |LT86j . Moreover, by the theory of linear elliptic equations with oblique boundary 
conditions (see jGTOlj ) and the assumed smoothness of data, one also has C 3,a (Q~) bounds for elliptic 
solutions, for any a < 1. 

6. Method of Continuity 

To complete the proof of Theorem 13.21 the standard method of continuity for nonlinear oblique boundary 
value problems is used, which is presented in | GT01I Section 17.9] and subsequently applied in | Urb97 | and 
| Urb01 |. This procedure was modified in | TW06 | in order to be applied to the Optimal Transportation Equation; 
and it is this method that is recalled here with trivial modification. Specifically, two procedures for applying the 
method of continuity were used in | TW06 |. The first method discussed there utilised foliations of both fi~ and 
Q+ to construct a family of suitable boundary value problems; this procedure will not be used here. The second 
method proved the existence of a function approximately satisfying the associated boundary condition, which 
subsequently allowed for the method of continuity to be applied without a domain variation construction. This 
is the procedure that will be recalled here (for the sake of completeness) with only trivial modification, as the 
Quotient Transportation Equation is constructed off the archetype of Optimal Transportation (see Subsection 

El- 

To start off, a key lemma from |TW06 | is recalled without proof: 

Lemma 6.1 ([TW06]). Let the domains Q~ and Q + and cost-function c satisfy the hypothesis of Theorem 
\3.1\ Then for any e > 0, there exists a uniformly c* -convex approximating domain, fi+, lying within distance 
e ofQ + , and satisfying the corresponding condition (|2.2ip for fixed 5q, together with a function u e C 4 (Q) 
satisfying the ellipticity condition (|2.8p , for f = a n ,i and the boundary condition (|l.lb|) for Q+. 

Remark. As /~(S n ) = r(a„j), the above lemma carries over trivially to the case of the Quotient Transportation 
Equation. 

With this in mind and making the denotation 

F H:=(f)^ [D 2 u-Dlc{-,T U )\, (6.1) 

the following family of boundary value problems are now defined: 

F[u t ] = tB{-,u t ) + (l-t)e lh - U0 F[u ], in (6.2) 

T Ut (Q") = Q+ (6.3) 

for t € [0, 1] with uo taken as the function indicated in Lemma 16.11 It is clear that u is the unique elliptic 
solution of (|6.2p - (|6.3p at t = 0. In this family of equations, is such that Qq = fi+ (corresponding to uq as 
depicted in Lemma [6~Tj) and = Q + . Given the assumed smoothness of c and the definition of T u depicted 
in (|2.17p . an e > can be chosen in Lemma l6~T1 small enough to guarantee that and are uniformly 
c-convex and c*-convex (respectively) relative to one another with corresponding uniform convexity constants 
independent of t, as fi~ and Q + are assumed to be uniformly c-convex and c*-convex (respectively) relative to 
one another. 

From Subsection l5.ll it is understood that the boundary condition (|6.3p is equivalent to the oblique condition 

G t ( ■ , Du) := 4>+(T u ( ■ )) = 0, on3fi^, 

where 4>f are defined for analogous to the construction of (p + for Q + in Subsection 15.11 From the obser- 
vations in the previous paragraph, it is clear that the family of boundary value problems correlate to uniformly 
oblique boundary value problems with a uniform constant of obliqueness independent of t. 
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To adapt the method of continuity from | GT011 Section 17.9], an a £ (0, 1) is fixed with Z set to denote the 
subset of [0, 1] for which t eT. implies the problem (|6.2p - (|6.3|) is solvable for an elliptic solution u t G C 2,a (Q7), 
with T Ut invertible. It is clear that the boundary condition (|6.3|) implies a uniform bound for Du t with respect 
to t. Upon noting that the inhomogeneity of ()6.2[) is uniformly bounded in t and satisfies (|2.22[) . f|2.23|) and 
(|2.24p , uniform estimates in C 2,1 (Q+) immediately follow, as all the solution estimates of Section [5] are clearly 
independent of t G [0, 1]. By compactness, it is then inferred that H is closed via the Heine-Borel Theorem. 
To show T. is open, the implicit function theorem is used along with the linear theory of oblique boundary value 
problems, as in | GT01I Chapter 17]. As H is open, closed and non-trivial, it is known that YL = [0, 1]; that is, 
there exists a unique elliptic solution u G C 3 (Q _ ) of the boundary value problem 

F[u] = B(-,u), in IT (6.4) 

T U (Q-) = Q+. (6.5) 

Thus, the part of Theorem 13.21 corresponding to unique solutions is proven. 



Remark. (|2.24|) guarantees uniqueness of the linearised boundary value problem via a straight-forward application 
of the Hopf boundary point lemma (see | GT01I Lemma 3.4]). Without the (|2.24p condition, one can not apply 
the method of continuity directly, as solutions will not be unique in general in this scenario. This will be discussed 
further in the next subsection. 

6.1. Application of the Leray-Schauder Fixed Point Theorem. To conclude this section, the Leray-Schauder 
fixed point theorem (see | GT01I Theorem 11.6]) will be used to relax the monotonicity criterion on the inho- 
mogeneity required by the method of continuity. The Leray-Schauder theorem has been used to similar ends in 
[LT861 . | Urb95| and | Urb01 |. 



Using the notation in (|6.ip , uq is defined to be the unique admissible solution of 

F[u ] = e°°, in QT , 
T U0 (Q-) = Q + . 

Theorem 13.21 indicates that such a unique uq exists. Moreover, by elliptic regularity theory, uq g C°°(Q~). For 
t G [0, 1] and ip G C 3 (Q~), the following family of problems are considered: 

F[u t ] = t(B(- , u Q + 1/0 + e u '-"°-^ - l) + (l - t)e Ut , in Or (6.6) 

T Ut (Q-) = Q + , (6.7) 

where B is only assumed to satisfy (|2.22p , (|2.23p and (|2.24l/vp . By Theorem 13.21 and elliptic regularity theory, 
for each i/j and t (|6.6p - (|6.7p has a unique admissible solution u t G C 3,a (Q~) for any a < 1. Consequently, the 
map T '■ C 3 (fi~) x [0, 1] — > C 3 (Q~) defined by T(ip, t) = u t — uq is continuous and compact; and T( • , 0) = 
for all ip G C 3 (fi~). If it can also be shown for t G [0, 1], that for all the fixed points of T( ■ , t) — that is, for 
any admissible solution u t of 

F[u t ] = tB{- , u t ) + (1 - t)e Ut , in rr (6.8) 
T Ut (Q-)=Q+ (6.9) 

— the estimate 

\\u t \\ CHTF) <C, (6.10) 

is satisfied with C independent of t, then by the Leray-Schauder fixed point theorem, T( • , 1) has a fixed point. 
Subsequently, this is equivalent to 

F[u] = B(-,u), infT, (6.11) 
r u (Q-) = Q+, (6.12) 

having an admissible solution u belonging to C 3 (Q~). As the solution estimates from Section [5] clearly apply to 
(|6.8p - (|6.9p independent of t, results from | LT86 | and linear theory for oblique elliptic boundary value problems 
can be applied to deduce (|6.10p . Subsequently, the monotonicity criterion on B may be relaxed to B z > 
relative to the previous existence criterion, thus finishing the proof of Theorem 13.21 
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Remarks. 

(1) As mentioned before, the admissible solution to (|6.11[) - (|6.12|) may not be unique as B z > is needed 
for the application of the Hopf Boundary point lemma. 

(2) It is recalled that the condition of B z > is used in the boundary C 2 estimate in Subsection 15.41 and 
thus, is still a required criterion, as estimates on solution of (|6.8p - (|6.9p need to hold independent of t. 

7. Conclusions 

While (|2.25ap - (|2.25b|) differ significantly from the Monge-Ampere equations coming from Optimal Trans- 
portation, the theory regarding their regularity is very similar. Specifically, the (A3w) represents both a necessary 
and sufficient condition for higher regularity of (|l.la|) - (|l.lb|) . The calculations in Subsection 15.31 indicate that 
the (A3w) condition is sufficient for Theorem 13.21 Given that the ellipticity criterion for (|l.lap - (|l.lb|) corre- 
sponds to solutions being c-convex, the counter-example that Loeper presented in [Loearj still applies in the 
current situation; and thus, (A3w) is also a necessary criterion for our higher regularity theory. 

Appendix A. c-convexity 

By using cost- functions, one can naturally generalise the classical notions of convexity. This generalisation 
is a powerful tool whose related notions and definitions will be used throughout the rest of this thesis. 

The following definitions have become standard in the literature surrounding Optimal Transportation (see 
|Caf96l [GM951 IMTW05I lvl\IQ8j for more detailed expositions on c-convexity). 

Definition 7.1. A c-support function of (p at xq is a function of the form c(x,y ) + a, where y G fi + and 
a = a(x , yo) is a constant, such that 

(p(x ) = c(x ,y ) + a, 

<p{x) < c(x,y ) + a, Vxerr. (7.1) 

Using the above definition, one may also define the notion of c*-support functions in a similar manner by 
switching x and y, and Now, the notion of support functions can be used to generalise the notion of 
convexity in the context of cost- functions. 

Definition 7.2. An upper semi-continuous function (p defined on is c-concave if for any point xo 6 Q~ , there 
exists a c-support function at Xq. Similarly, an upper semi-continuous function ip defined on fi + is c*-concave 
if for any point y £ Q + , there exists a c* -support function at yo. 

From these definitions, one can also easily ascertain the notion of a c-convex function by simply switching 
the direction of the inequality shown in (|7.ip . The definition for a c*-convex function is also thus obtained 
by replacing c-support functions with c*-support functions in Definition ^. 2l and reversing the inequality in (|7.ip . 

Before proceeding onto other definitions, some properties of c-concave functions can be readily observed 
through facts in classical analysis. As it is assumed that the cost-function c is smooth, any c-concave function 
(p is semi-concave; that is, there exists a constant C such that (p(x) — C\x\ 2 is concave. It is a classical result 
that now shows <p is twice differentiable almost everywhere. In addition to this, it is readily demonstrated (via 
a Perron process) that if (cpk) is a sequence of c-concave functions and (pk — > (p, then <p is c-concave. Clearly, 
there are analogous results for c-convex and c*-convex/concave functions. 

Remark. It is important to note that the above definitions correspond with the classical notions of concavity in 
the special case where c(x, y) = x ■ y. In this situation, a c-support function is simply a support hyperplane at 
a particular point. This special case provides a good basis of intuition for more general cost-functions. 

In addition to using generalised definitions for concavity of functions, one also needs to consider analogous 
generalisations of concavity with respect to sets. In particular, such notions will be of key importance in 
the construction of barriers for boundary gradient estimates, in addition to being necessary for making the 
statements of current regularity theorems in Optimal Transportation. With that, the notion of c-segments 
(which generalise the notion of line-segments in classical convex analysis) is next defined: 
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Definition 7.3. A c-segment in Or with respect to a point y is a solution set {x} to D y c(x,y) = z for z on 
a line segment in W. A c*-segment in Q + with respect to a point x is a solution set {y} to D x c(x,y) = z for 
z on a line segment in R" 

By conditions (Al) and (A2), it is clear that a c-segment is a smooth curve; and for any two points x , x x G R" 
and any y G R", there exists an unique c-segment connecting xo and Xi relative to y. 

Definition 7.4. A set E~ is c-convex relative to a set E + if for any two points x 0l x% G E~ and any y G E + , 

the c-segment relative to y connecting xo and Xi lie in E + . Analogously it is said that E + is c*-convex relative 
to E~ if for any two points yo.yi G E + and any x G E~ , the c* -segment relative to x connecting yo and yi lies 
in E+. 

Remarks. 

(1) In general the notion of c-convexity, with respect to sets, is stronger than the classical notion of 
convexity. This is most readily observed when taking a ball for the set in question. Indeed, a ball may 
not be c-convex (relative to another given domain) at an arbitrary location. However, a sufficiently 
small ball will be c-convex if c is C 3 smooth. In light of this example, it is insightful to consider an 
alternate statement of Definition 17.41 

E~ is c-convex with respect to E + if for each y G E + , the image D y c(E~ , y) is convex in 
R". 

(2) One may reconcile the above alternative definition with (|2.2ip by pulling back that equation by the 
mappings c y (-,y) to see that this immediately implies that the image domains of c y (Q~,y) are 
uniformly convex with respect toy £ Q + | TW06| . 

(3) An analogous statement may be made for a set that is c*-convex relative to another set by considering 
c*(x,y) = c(y,x). 
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